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Abstract 

A magnetized large-volume liquid scintillator might be used as a tracking 
detector with charge identification to measure high-energy neutrinos. I con- 
sider neutrino oscillations for a beam from a low-energy neutrino factory to 
such a detector for the baseline CERN-Pyhasalmi, with the length of 2288 km. 
The range to study the oscillation parameters extends to sin^ 29i3 ^ 10~^ for 
reasonable assumptions for the detector and the beam. 

1 Introduction 

I study the low-energy neutrino factory [H El [3l IH O El [7] with the 2288 km long 
baseline CERN-Pyhasalmi, using a magnetized large volume liquid scintillator. As 
the prototype model I take a magnetized LENAP [9l[l0l[IIl[ia[l3l[H[ISlIl6],a 
30 m wide cylinder in vertical orientation, by default with 50 kton fiducial mass. 

It has been suggested that a large volume liquid scintillator indeed has capacity 
to measure high-energy neutrinos |T71 [l8j. On-going studies hint that with a good 
photodetection a magnetic liquid scintillator can be used to measure the charge of 
muons. The performance may depend on the design of the detector and on the 
choices of the components. Here I assume an optimistic performance and study 
different options, to set the goals for the design of the experiment. 

Previously a conventional beam of 1-6 GeV from CERN to LENA at Pyhasalmi 
mine has been considered[19j, as well as high-energy beta beams [20]. The conven- 
tional beam is a very viable option to study neutrino parameters if the third mixing 
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Figure 1: Spectrum of neutrinos from a 5 GeV muon beam. 

angle is not too small, i.e. sin^ 26^13 > 3 • 10~^, and the beta beam with the same 
setup might extend the range to sin^ 26^13 > 1 • 10"^. With a larger detector or 
stronger beams those limits might be extended by a factor of 2-5. 

The Pyhasalmi Mine provides a 1400 m (4000 m.w.e.) deep underground site in 
the middle of Finland. The density profile for the CERN-Pyhasalmi baseline has 
been modelled well [2J4 |22j, and the accuracy of the average density can be taken 1 
%. In these simulations, a simplified profile is used [23j . 

2 Low-energy neutrino factory 

For the CERN-Pyhasalmi baseline the first oscillation peak is at about 4.2 GeV with 
usual parameters. Hence a beam with an energy range 1-5 GeV appears the most 
interesting for this baseline. 

Here I study neutrinos made in a neutrino factory by the decay of high-energy 
muons: /x~ z/^e~I7g, and respectively for the positive charge. The golden signal is 
the wrong-sign muon appearance, as well as the respective disappearance channels. 
A wrong-sign electron appearance would be more challenging (platinum channel). 
The tau appearance (silver channel) is relevant for higher energies. 

For the detector a technology to measure the neutrino helicity is required. In 
practise the lepton charge is measured using magnetic field. Non-magnetized detec- 
tors with a low-energy neutrino factory were studied in [24] . 

The technology for the neutrino factory is currently under development. It is 
not yet completely clear how feasible it will be, but in any case the low energy 
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neutrino factory (5 GeV) is expected to be substantially cheaper than the ultimate 
high-energy neutrino factory (20-50 GeV). 

I assume for the beam power 1.1 • 10^^ usable muon decays per year, and 5 year 
running time for both polarities. 

The detector capable to measure neutrinos from the neutrino factory will also 
work well for conventional or beta beams. While the beta beam may be rather 
a competitor for the neutrino factory, a conventional beam could be considered 
complimentary, or rather a predecessor to the neutrino factory, given that the cost 
of the beamline will be a tiny fraction the cost of the neutrino factory. Probably 
the neutrino factory would be built only after the preceding beam fails to measure 
the interesting quantities due to 9is being out of range. Hence I assume here that 
there will be an accompanying conventional wide band beam, as described in [19j . 



3 Detector 

As the detector model I take a magnetized LENA. The plain LENA is planned to 
consist of 50 kton of liquid scintillator, in a vertical cylindrical tank with the radius 
of 12 m and the height of 100 m. It is surrounded by a buffer of 2 m, and outside it 
is shielded by water (totally 100 kton). The scintillation light is read by more than 
10 000 photosensors. 

I assume that the detector can be magnetized, either surrounding the tank with 
a circular coil to generate a homogeneous vertical field inside the tank (and outside a 
inhomogeneous cylidnrically symmetric field in the opposite direction) or otherwise. 

With a good photodetection system the large volume liquid scintillator can be 
used as a tracking detector. One or more tracks can be reconstructed from the 
measured photon time profiles. With a good reconstruction one can define the lepton 
fiavor and estimate the neutrino energy. This requires a fast scintillator (3 ns) and a 
good time resolution for the photosensors (better than 3 ns), as well as an ability to 
count multiple photons arriving in a single photosensor with a reasonable accuracy. 
While the considered concept is formally a totally active scintillator detector, it is 
fundamentally diflFerent from the segmented scintillator [6j which the term TASD 
usually refers to. 

The lepton charge can be defined if its curvature under magnetic field is mea- 
surable. Initial studies point that the required magnetic field to cause an observable 
bending is about 0.02-0.05 T, depending on the time properties of the scintillator 



and the photosensors, and also on the event topology (See Fig. 10). On the other 
hand, the magnetic field must be strong enough to separate the magnetic curvature 
from stochastic deviations, mostly multiple Coulombian scattering. The scale of the 
multiple scattering is defined by the scattering length, in liquid scintillator 42 cm, 
compared with 1.8 cm in iron. For typical track lengths the magnetic bending is 
equal to the mean stochastic deviation with B = 0.02 T, and to get suflBcient sta- 
tistical confidence at least 5 = 0.1 T is required. Further simulations are required 
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to define the optimal magnetic field to reach a sufficiently low misidentification 
probability, particularly because the multiple scattering is not Gaussian. 

The magnetization requires photosensors that work under such a magnetic field. 
The availability of such detectors of sufficient time resolution and dynamic range is 
not yet clear, but we have several years to develop such devices before the decisions 
for the neutrino factory are due. Decisions to construct LENA may be, however, 
due earlier than that, and for that the selection of the photosensor technology is a 
rather urgent task, as well as the decision to magnetize the detector (or at least to 
keep an option for it). 

Other than the magnetization, there are some chances to distinguish the neutrino 
helicity by other methods, though only statisticallv[24j. The most significant may be 
identifying the recoil nucleon — the identification is rather clear but due to possible 
intranuclear interactions the nucleon flavor may be changed at 10-20 % probability, 
and the presense of charged pions with unidentified charge spoils the possibility 
altogether. Given the relatively moderate cost of the magnetization compared with 
the total cost of the experiment, using a non-magnetized scintillator may not be a 
cost-effective way to go, unless the magnetization turns out to be unfeasible due to 
unavailability of suitable photosensors or other reasons. For the platinum channel 
(electron appearance), however, the given statistical methods may be the only way 
to get charge information unless the photodetection can be improved beyond present 
assumptions (e.g. by photographic methods). 

I consider also a larger size, 150 kton, for comparison. This could be achieved by 
enlarging the size, using the outer shield as additional volume or building another 
dedicated beam detector, probably aligned towards the beam. Increasing the beam 
power (which is a very unknown quantity) would have a similar effect as increasing 
the detector size. Given that the large-volume liquid scintillator is expected to 
be much cheaper than a segmented scintillator or iron calorimeter (and a dedicated 
high-energy detector costs only a fraction of an all-range detector), in this case larger 
detectors are more affordable and even desirable to balance the total investment. 
Anyway, the detector will have also a research programme of its own, for atmospheric 
neutrinos, proton decay, cosmic rays and low-energy neutrinos. 

Throughout this work I assume the energy resolution of 5 %, unless otherwise 
stated. The resolution depends quite substantially on the properties of the scintilla- 
tor, phototubes and electronics, but according to on-going work, this assumption is 
realistic, unless the photoaccuracy is spoiled by the magnetisation. From the anal- 
ysis point, the magnetization itself just improves the accuracy of energy definition. 

The detection aspects will be described in more detail elsewhere. 

4 Simulations 

The simulations are done using the GLOBES codebase [25l [26] . embedded within an 
own code. The AEDL file described in [27j was used as a starting point. 
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Figure 2: The CP range with a combined set of a neutrino factory and a wide band 
beam for the same 2288 km baseline. Left 50 kton and right 150 kton. There is a major 
difference, the larger is much better. The colors represent x values, so that yellow is 
about 3(7, red 2a and blue la. 



The hottest motivation for the long basehne oscillation experiments is to discover 
leptonic CP violation. This is expressed with the parameter 5, values differing from 
or TT meaning CP violation. The observable range of the CP angle 5 for the 
assumed beam and baseline is shown in Fig. [2} with both 50 kton and 150 kton 
fiducial masses. At best the range extends to sin^ 2^13 ~ 8 • 10"^ with 150 kton, 
with 50 kton sin^ 2^13 ^ 2 • 10~^. This is significantly better than that achievable 
with conventional beams or beta beams, though a fair comparison depends on the 
achievable beam powers. 

The considered baseline is particularly suitable to study the mass hierarchy, i.e. 
the sign of 611223. The neutrino factory extends the search of the mass hierarchy 
to smaller mixing angles than any other beam. In the given examples, with 150 
kton detector sin^ 2^13 ^ 3 • 10""^ is reachable for the worst 5, and for the best case 
sin^ 2^13 ^ 8 • 10~^. Without the assumed conventional beam the performance curve 
at 5 plane has glitches towards lower miging angles, i.e. the conventional beam 
improves the performance at lower mixing angles. 

Measuring the last unknown mixing angle, 9 is is also one of the main goals 
of the experiment. If other experiments measure it before neutrino factory, the 
neutrino factory will be rather a precision machine than discovery machine (though 
measuring a large mixing might detract seriously from the motivation to build the 
neutrino factory). The range to measure 6^13 extends to sin^ 26^13 ^ (2 . . . 9) • 10"^ 
with 50 kton, or sin^ 2^13 ^ (0.8 . . . 3) • 10~^ with 150 kton. Without the conventional 
beam the reach with the worst 6 is poorer. 

To set the goals for the detector design I studied the performance for diflFerent 
detector properties. In figure [6] is shown as example the capacity to measure the 
mass hierarchy with diflFerent relative energy resolutions. That shows that a 20 
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plotfile using 1:2:3 
required fiducial nass of the detector to define nass hierarchy 



plotfile using 1:2:3 
required fiducial nass of the detector to neasure non-zero theta_13 
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Figure 3: The capacity to measure the mass hierarchy (left) and the 9is reach (right), 
with both a neutrino factory and a wide band beam together. In the horizontal axis is 
the fiducial mass of the detector, used as a scaling factor including also variations in the 
beam power, the detector efficiency and the running time. Here and in the other similar 
figures S = 37r/2. 
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Figure 4: The capacity to measure the mass hierarchy in the plane of 6^13 and S, with 
50 kton and 150 kton detectors and neutrino factory and conventional beam together. 
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theta reach; theta_13 and delta for 5B kt 



plotfile using l;2;3 
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Figure 5: The reach of 6^13 for a 50 kton detector and for a 150 kton detector with a 
neutrino factory and a wide band beam together for the same baseline. 



% resolution is required in order not to spoil the capacity to measure the mass 
hierarchy, and 5 % resolution is sufficient. 
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Figure 6: The capacity to measure the mass hierarchy for shown relative energy resolu- 
tions (150 kton detector for neutrino factory only), left 5 = 0.3 and right 5 = 2.2. The 
assumed 5 % resolution is rather well-suited, as better resolution will not improve the 
capacity very much. 



I also studied the inclusion of the platinum channel with statistical cahrge recog- 
nition. As an example the results for CP violation are shown in Fig. (Q) Runs with 
50 % (pessimistic) 70 % (realistic), 90 % (optimistic) and 99 % (utopic) charge 
identification probabilities were made, showing an improved performance at larger 
9is with a good electron charge recognition. A plain statistical charge identification 
is of little benefit. 
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5 Conclusions 



Evidently a magnetized large volume liquid scintillator can be used as a far detector 
for the neutrino factory. For a low-energy (max. 5 GeV) beam a magnetized LENA 
with magnetic-field-tolerant photosensors would be suitable, higher energies need 
a dedicated detector aligned in the beam direction. While the 50 kton detector 
performes well, the vertical orientation being a small burden, a larger detector would 
be more optimal. The large volume scintillator is probably cheaper than a segmented 
scintillator and much more accurate than an iron calorimeter and has a rich research 
programme of its own in particle astronomy and nucleon decay searches. 

To reach the required accuracy an appropriate photodetection system working 
under magnetic field must be designed. Max a few-ns time resolution and multi- 
photon capacity are needed. Additional studies are needed to prove whether this is 
realistic. Further simulations are also required to settle the optimum magnetic field, 
initial studies suggesting that 0.1 T may be enough. 

The low-energy neutrino factory with a 2288 km baseline was found to have 
a good capacity to study the neutrino properties. Neutrino oscillations up to 
sin^ 2^13 ^ 10"^ . . . 10"^ can be studied. For the lowest 6^13 this beam outperforms 
the beta beam and any superbeam for a similar baseline. For a fair comparison, 
however, one has to take into account that at the cost of the neutrino factory one 
could build a whole lot of detectors for a conventional beam. 
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Figure 7: For comparison, the reach in 9is with only a neutrino factory, 50 and 150 kton 
detectors. 



paraneter area to discover CP violation 
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Figure 8: For comparison, the CP range with only a neutrino factory with a 50 kton 
detector. Left for the horizontal layout (no high-energy cut) and right for a vertical 
layout, i.e. reduced high-energy efficiency. The difference between these is barely visible, 
but the shapes differ remarkably from Fig. [2] that includes the conventional beam, too. 
For small 9is and small CP violation the adjoint conventional beam is very important. 



10 





Figure 9: The capacity to study the CP violation with the neutrino factory only, including 
the golden and the platinum channel, with 150 kton. The figures are for 50 %, 70 %, 90 
% and 99 % charge resolution (the first is no resolution at all). Evidently the statistical 
charge resolution brings little improvement. 
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Control Options Event Measurement Analyse View Layer Help 




COMMAND: Check charge 
event generated 
^/lEW: top 
LAYER: photons 



fRSLVERDICT 
Error in measured energy 4.05 MeV= 0.1 4 % 
Error in lepton energy-21.33 MeV= 0.80 % 
Error in lepton track 0.09 m = -1 %, vertex: 0.1 m. 
Error angle L 0.08 rad = 5 deg (p 0.43 rad = 25 deg) 



DETECTOR Volume = 21 206 m3 B = 0.05T 

Photosensor coverage = 6 % PDE of photosens brs = 1 0O % I 
ORIGINAL EVENT QE with neuthno energy 3000.0 MeV 
Depositable energy 2858.60 MeV Measurable energy 2963.00 MeV 
muon:2664.1 5 MeV and 1 2.S7 m. vertexEnergy =0.00 MeV 
proton:1 94.45 MeV and 0.290 m. 



MEASUREMENT 

measured 307529 photons of 31 .1 6 M. (0.99 %) 
FIT (studied charge) 
ln(L) = 1109822 s=0.00 

Vertex at (-3. 41, 0.02, 0.1 3)0.00 MeVt0 = 60.79 ns. 
Deposited energy 2862.65 MeV Measured energy 2967.05 MeV 
Inferred neutrino energy 31 01 .45 MeV with uncertainty 1 34.39 MeV 
Neutrino energy from lepton angle: 2949.58 MeV [QES] 
[0] muon(-1)2643 MeV and 1 2.78 m. 
[1] proton:220 MeV and 0.360 m. 
Predict 307745 photons of 2.28 M emitted. (1 3.51 %) 
best fit original, with measured E= 2963.00, Chi = 922337203685477580 




- I PMT quality perfect 

1 77 photons 1 71 photons 

Duration 82.94 ns Duration 60.00 ns 

Time 60.88 ns Time 61 .13 ns 

Mean time 6.41 ns Mean time 6.1 6 ns 

Decay signal = 1 Decay signal 0.76 

+1 +0.76 

ln(L) 0.14 



Figure 10: A screenshot from a simplified simulation program, with B = 0.05 T, showing 
how to resolve the curved track for a multi-GeV muon using the photon observation time 
information of all individual photosensors, with an optimal photodetection system. The 
blue line is the original muon track and the red one the fit. 
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